On the ground state of a completely filled lowest Landau level in two dimensions 
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There exists a widely believed opinion, that the many-body ground state of a two-dimensional 
electron system at a completely filled lowest Landau level (the filling factor v = 1) is described by 
the so-called Hartree-Fock wave function, and that this solution is the unique, exact eigenstate of 
the system at v = 1. I show that this opinion is erroneous, construct an infinite number of other 
variational many-body wave functions, and discuss the properties of a few states which have the 
^ , energy substantially lower than the energy of the Hartree-Fock state. 
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The nature of the ground state of a system of two-dimensional (2D) interacting electrons in strong magnetic fields 
B was a subject of intensive investigations in the past years p|-p4[. A great deal of attention was given to a partially 
filled lowest (y < 1) fli]-[T(| or higher Landau levels (v ^ 1) |17| |2^. The case of a completely filled lowest Landau 
level v = 1 was not adequately investigated in the literature. The only trial many-body wave function for this v, the 



so-called Hartree-Fock solution, 



*™ =-^det|^(r i )|, L 3 =0,1,...,A-1, (1) 



e HF = -tt/2 = -1.57080 (3) 



was proposed in Ref. p4|. This is a Slater determinant, constructed from the lowest-Landau- level single-particle states 

V>4r) = ^2i=exp(-zz72), (2) 

o . 

and it is assumed that in the A^-electron system the states L = to L = N— 1 are occupied by electrons. Here 
z = (x + iy)/X is a complex coordinate of an electron, A 2 = 2l 2 = 2hc/eB, I is the magnetic length, L — 0, 1, 2, ... is 
| the angular momentum quantum number, and h, c, and e are the Planck constant, velocity of light and the electron 
charge, respectively. The wave function (Q) is an eigenfunction of the kinetic energy operator, with the eigenenergy 
E = NTiuj c /2. It coincides with the Laughlin many-body wave function || at v = 1. The state (|l|) is characterized, 
in the thermodynamic limit, by a uniform 2D electron density n s (r) = 1/irX 2 and the energy per particle j(| 
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| in the ^-independent energy units e\ 

There exists an opinion, that the Hartree-Fock many body wave function (Q) is the only one possible solution of 
the problem at v = 1. Although this opinion is incorrect, it seems to be widely believed. The aim of this Letter is to 
show that the Hartree-Fock many-body wave function is not the only one possible trial wave function for the ground 
state of a 2DES at v = 1, and to demonstrate a number of other variational solutions of the many-body Schrodinger 
equation at v = 1, which have the energy substantially lower than the energy of the Hartree-Fock state. 
q I argue in three different manners. 

The first argument is quite simple: in order to prove that the function (Q) is not unique at v = 1, it is sufficient to 
present another explicit example of a trial wave function. Consider for instance the Wigner crystal many-body wave 
K~2 ' f unct i on @7 
■ 

= det \ XL= o(Ti, R,-)l, (4) 



XiO-i.Rj) =M^-~R J )e-^< B * R ^\ (5) 

Here (f>o is the flux quantum and Rj are points of a triangular lattice [Q, distributed over the 2D plane with the 
average density n s . This function depends on magnetic field and can be considered at v = 1. One can easily see 
that the projection of the Wigner crystal wave function onto the Hartree-Fock one is not unity. Expand the function 
Xo( r , R) in a set of the lowest-Landau-level eigenstates ipL(r), 
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Xo(r,R) = 
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where C L (Z) = Z L exp(-ZZ* /2)/VlA. Then, for the two-electron function ^ I get 
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All the basis functions det \ipLj ( r i)|/v2! here are orthonormal. One of them (with L\ = and L 2 = 1) is the 
Hartree-Fock function (nh. The projection 
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(8) 



is evidently smaller than one. The same derivation can be easily performed at any N . In the thermodynamic limit 
N — > oo the projection (||) tends to zero, P — * 0. 

Arbitrarily varying the vectors Rj (square lattice, other types of the lattice, random distribution), as well as taking 
the functions 



i 



det |xi(ri,Rj-)|, 



(9) 



with other angular momentum index L, one can easily get an infinite number of other explicit examples of many-body 
wave functions, different from (Q). 

The second argument is based on a standard degenerate perturbation theory. Consider the problem of N 2D 
electrons in a perpendicular magnetic field, in the presence of a neutralizing positive background, which has the form 
of a disk with the radius R and the charge density +en s . The electroneutrality condition requires that 7rP 2 n s = N. 
Assume that the total Coulomb energy of the system (electron-electron plus electron-background plus background- 
background interaction energy) can be considered as a perturbation. The ground state of the unperturbed problem 
is highly degenerate, therefore, one should use a degenerate perturbation theory. In the one-electron problem, one 
should search for a solution in the form v]/[ Ar = 1 ] (r) = 5~J£L CltPl{y)- This expansion contains all L-terms, from L = 
to L — oo. In the A^-electron problem, one should write, similarly, 



det|V>L,(r;)|, 



(10) 



where, again, all Li vary from Li = to Li — oo. The Hartree-Fock Slater determinant is the only one term in this 
expansion (L\ = 0, Z>2 = 1, • • Ln = N — 1). Obviously, any variational solution of the many-body Schrodinger 
equation can be also searched for in the form of an arbitrary linear combination ([h]). Hence, the number of all 
possible trial wave functions (at any v) is not one but infinite. 

Finally (the third argument) I discuss the problem in the finite-cell geometry. In a number of papers (see e.g. 
Ref. the problem of a single 2D electron in a perpendicular magnetic field was considered in a finite rectangular 
cell (0 < x < a, < y < b), with periodic boundary conditions at the boundaries of the cell. According to Q, the 
"boundary condition requires that ab/2-!rl 2 be an integer m", and "there are m different single-electron states in the 
cell" 
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where 1 < j < m and Xj = 2irl 2 j/b. From this statement it follows that, at v = 1 the number of electrons per cell 
is exactly equal to the number of single-particle states, and therefore there exists only one possible way to construct 
the many-body wave function. Consider this argument in some more detail. 
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Consider a rigorous mathematical formulation of the single-electron problem in a finite-cell geometry. The wave 
function sought should satisfy the differential equation 



-h 2 d 2 x + (-ihdy + eBx/cf <f>(x,y) = 2mE(j)(x,y), (12) 
in the area < x < a, < y < b, and the boundary conditions 

<fi(x, y = 0) = cj)(x, y = b), d y 4>{x, y = 0) = d y (j)(x, y = b), at all < x < a, (13) 

<f)(x = 0,y)= 4>(x = a, y), d x (j)(x = 0, y) = d x (j)(x = a, y), at all < y < b. (14) 



Standard substitutions, <j)(x,y) = (p(x) exp(ik y y), £ = x/l + k y l, e = 2E/Tiuj c , lead to the following equation for 
(p{£) = ip{x = l£,-k v l 2 ), 

The second-order differential equation ( |l5| ) has two independent solutions, e.g. 

e) - e-« 2 / 2 [1 + (1 - e)e/2l + (1 - e)(5 - e)^/^ + ■ ■ ■} (16) 

and 

*a(e, e) = e"« 2 / 2 [£ + (3 - e)£ 3 /3! + (3 - e)(7 - e)£ 5 /5! + ...]• (17) 
The total solution of eq. ( |l2|) is then written in the form 

4>{x, y) = e lk ^ [C'Mx/l + k y l, e) + C 2 <f> 2 {x/l + k y l, e)] , (18) 



with two arbitrary constants C\ and Ci. The boundary condition (|13| ) requires that k y — 2irm/b with integer to. If 
the second boundary condition was imposed at infinity [4>{x — ±oo,y) — 0] one would get a conventional solution of 
the problem, with e n = 2n + 1, n = 0, 1, . . ., and Landau eigenfunctions. In the finite-cell geometry, the boundary 



conditions (14) require that 

{$>i{k y l, e) - $i(a/I + k y l, e)}d + {3> 2 (k y l, e) - $ 2 (a/^ + e)}C 2 = 0, . . 

Equations ( |Tgj ) determine the spectrum of eigenenergies e = e n (a/l,b/l), n — 0,1,..., and eigenfunctions of the 
problem (they are different from the Landau ones). Obviously, the boundary conditions do not restrict the number 
of eigenstates, and the single-particle spectrum remains unlimited. 

One can easily verify, by means of a direct substitution, that the wave functions (|ll|), proposed in do not satisfy 
the boundary condition ([l4|). Instead, they satisfy the boundary condition (f>(x — a,y) — exp(iay /l 2 )(f>(x = 0,y)- 



Hence, the functions ( |ll[ ) are not the eigenfunctions of the boundary- value problem (|12|)-(|14|), and the conclusion of 
Ref. |g| about a finite number of possible single-particle states in the considered problem is incorrect. 

Now I present results of calculations of the energy of a few trial wave functions, different from the Hartree-Fock 
solution ([j]). I consider the functions eq. ([}]), with L = to 5, and with a triangular lattice of points Rj, uniformly 
distributed over the 2D plane with the average density n s . The states ^ l describe the properties of the system at all 
Landau-level filling factors v < 1, continuously as a function of magnetic field. Like in Pq| , I calculate the Hartree 
energy exactly, in the thermodynamic limit, and the exchange-correlation energy approximately, at a finite number N 
of lattice points. Thus calculated energy per particle is the upper limit to the true energy cl, £l < ^l(N) [ p5| . 

Figure [j] exhibits the energy cl(N) of the states ^l, at v = 1, as a function of the number of lattice points 
N involved in calculation of the exchange-correlation energy. All the states ^l, L = 0, 1, . . . , 5, have the energy 
lower than the Hartree-Fock one. Some of them (L = 1, 3, 4 and 5) have the energy lower than the classical Wigner 
crystal |2{J. The energy of the state ^^=3 estimated with N = 187 lattice points, is more than 46% lower than the 
energy of the Hartree-Fock state. All the considered states are characterized by a strong overlap of the neighbor 
single-particle wave functions and by a uniform 2D electron density n s (r) = 1 /ir\ 2 at v = 1 (see for example, Figure 3 
from Ref. ]25|). Subsequent investigations of other trial wave functions (with other L and/or different symmetries of 
the lattice Rj) may lead to even lower estimates of the ground state energy of the 2DES at v = 1. A huge variational 
freedom of the functions *$>l, which has been clearly demonstrated in this Letter, opens up wide possibilities to search 
for better approximations to the ground state, both at v — 1 and at all other v. 
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FIG. 1. Calculated energies el(N) per particle of the states (in units e 2 y / n7), as a function of the number of (triangular) 
lattice points N, for the states L — to 5 at a completely filled lowest Landau level v = 1. The Hartree contribution is 
calculated in the thermodynamic limit, the exchange-correlation contribution - for a finite number N of lattice points. For 
comparison, the energies of the Hartree-Fock state (the thermodynamic limit, TDL) and of the classical Wigner crystal |2(]] 
are shown by solid and dashed lines respectively. 

This work was supported by the Graduiertenkolleg Komplcxitat in Festkorpern, University of Regensburg, Germany 
and the DFG-Sonderforschungsbereich 348 (Nanomctcr-Halblcitcrbauclemente). I thank Ulrich Rossler, Nadejda 
Savostianova, Vladimir Volkov, Vladimir Sandomirsky, Oleg Pankratov, Robert Laughlin, Pawel Hawrylak and Allan 
MacDonald for discussions and useful comments. 



Fukuyama H., Platzman P. M. and Anderson P. W., Phys. Rev. B, 19 (1979) 5211. 
Yoshioka D. and Fukuyama H., J. Phys. Soc. Japan, 47 (1979) 394. 
Yoshioka D. and Lee P. A., Phys. Rev. B, 27 (1983) 4986. 
Maki K. and Zotos X., Phys. Rev. B, 28 (1983) 4349. 

Yoshioka D., Halperin B. I. and Lee P. A., Phys. Rev. Lett., 50 (1983) 1219. 

Laughlin R. B., Phys. Rev. Lett, 50 (1983) 1395. 

Haldane F. D. M., Phys. Rev. Lett, 51 (1983) 605. 

Yoshioka D., Phys. Rev. B, 29 (1984) 6833. 

Lam P. K. and Girvin S. M., Phys. Rev. B, 30 (1984) 473. 

Levesque D., Weis J. J. and MacDonald A. H., Phys. Rev. B, 30 (1984) 1056. 

Haldane F. D. M. and Rezayi E. H., Phys. Rev. Lett, 54 (1985) 237. 

MacDonald A. H., Aers G. C. and Dharma-wardana M. W. C, Phys. Rev. B, 31 (1985) 5529. 
Kivelson S., Kallin C, Arovas D. P. and Schrieffer J. R., Phys. Rev. Lett, 56 (1986) 873. 
Morf R. and Halperin B. I., Phys. Rev. B, 33 (1986) 2221. 
Jain J. K., Phys. Rev. Lett, 63 (1989) 199. 

Halperin B. I., Lee P. A. and Read N., Phys. Rev. B, 47 (1993) 7312. 

Aleiner I. L. and Glazman L. I., Phys. Rev. B, 52 (1995) 11296. 

Koulakov A. A., Fogler M. M. and Shklovskii B. I., Phys. Rev. Lett, 76 (1996) 499. 

Fogler M. M., Koulakov A. A. and Shklovskii B. I., Phys. Rev. B, 54 (1996) 1853. 

Moessner R. and Chalker J. T., Phys. Rev. B, 54 (1996) 5006. 

Fogler M. M. and Koulakov A. A., Phys. Rev. B, 55 (1997) 9326. 

Jungwirth T., MacDonald A. H., Smrcka L. and Girvin S. M. , Physica E, 6 (2000) 43. 
MacDonald A. H. and Fisher M. P. A., Phys. Rev. B, 61 (2000) 5724. 
Bychkov Y. A., Iordanskii S. V. and Eliashberg G. M., JETP Lett, 33 (1981) 143. 



Mikhailov S. A., A new approach to the ground state of quantum- Hall systems, preprint cond-mat/0003228 
Bonsall L. and Maradudin A. A., Phys. Rev. B, 15 (1977) 1959. 



4 



